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Abstract. - We present a simple model of electrical transport through a metal-molecule- 
metal nanojunction that includes charging effects as well as aspects of the electronic structure 
of the molecule. The interplay of a large charging energy and an asymmetry of the metal- 
molecule coupling can lead to various effects in non-linear electrical transport. In particular, 
strong negative differential conductance is observed under certain conditions. 



Introduction. Single molecule electron transistors offer exciting perspectives for further 
minituarisation of electronic devices with a potentially large impact in applications. To date 
several experiments have shovirn the possibility to attach individual molecules to leads and to 
measure the electrical transport. Two terminal transport through a single molecule or 
other nanoscale objects has been achieved by deposition of the object between two fixed 
electrodes or a conducting-tip STM above an object attached to a conducting substrate 
Interesting and novel effects, such as negative differential conductance (NDC), were observed 
in one experiment (l2| , which still needs satisfactory theoretical explanation. 

Several factors are important for single-molecule transport: For nanoscale objects the ca- 
pacitance C is very small. Consequently, the energy to charge (or uncharge) the molecule 
Ec = e^/2C can be very large, of the order of electron volts. This leads to the phenomenon 
of Coulomb blockade and makes room temperature single-electron transistors (SET) based on 
such molecules a distinct possibility. In contrast to SETs based on metallic islands [p^p^, 
molecular devices have a more complicated electronic structure that, in principle, can be chem- 
ically 'designed'. Gaps in the I-V curve are not only determined by Ec but predominantly 
by the structure of the molecular bands Therefore, it is important to consider the in- 

terplay of charging effects with the specific structure of the molecular orbitals. For electronic 
transport we will see that, in particular, the specific spatial structure of the molecular orbitals 
is crucial. 

In this letter we study the impact of such a low-energy electronic structure on electronic 
transport and demonstrate that it can result in non-trivial conductance under reasonable 
generic assumptions. A full quantitative treatment of a molecule in contact with metal elec- 
trodes including many-body interactions on the molecule and strong molecule-electrode cou- 
pling is still out of reach. In this paper, we therefore study a simple model for a 'generic' 
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molecule that can be derived in principle from full electronic structure calculations and that 
includes capacitive interactions within the Coulomb blockade model as well as absorption and 
emission of photons. In contrast to most other theoretical work 
coupling between electrodes and the functional part of the molecule, a scenario that can be 
realized by proper design of molecule and contacts. We predict that the interplay of charg- 
ing effects, the spatially non-trivial electronic structure and resulting asymmetric coupling to 
the electrodes leads to interesting phenomena in non-linear transport. We find that the I-V 
curve can show strong NDC behaviour and thus may provide a possible mechanism to explain 
recent experiments p2| showing similar NDC behaviour. Finally, we give a specific example 
that displays the required electronic structure in a simple aromatic molecule. 

The Model. In the low-voltage regime only a few of the molecular levels will contribute to 
transport. For the purpose of this paper, we assume without loss of generality that there are 
only two participating molecular levels that are both unoccupied at zero voltage, which we 
designate as LUMO and LUMO+1. In particular in tt electron systems of aromatic molecules 
it is relatively easy to realize a situation, where two closely spaced MOs are separated far 
from both from the HOMO and the other LUMOs. The other MOs are assumed inert (always 
occupied or always empty), as we ignore co-tunnelling effects. 

The MO Hamiltonian of the 'reduced' system can be written as: 

n = ^ e,^c|^C,^ -1- V^jklc\^c]^,Cka'Cla, (1) 

i(T ijklaa' 

where the operators Cio-,c|^ destroy /create electrons with spin a in MO i. This Hamiltonian 
contains a large number of parameters even for the two level system. In order to make contact 
with the quantum dot physics and to elucidate similarities and crucial differences that occur 
in single molecule transport it is useful to introduce an effective model for the molecular 
Hamiltonian: 

H^oi = eiN, + e2N2 + Ec {N, + N^f (2) 

+^Y.^i{N, - 1) + ^ 5:^'iA'Ci.-2., 

/ (T.{T' 

where Ni is the occupation of the MO I. As unit of energy we take Ae = €2 — ei, the "bare" MO 
splitting. The diagonal electronic repulsion terms of Eq. (]l|) for the two MOs together with 
capacitive interactions with the leads can be rewritten as Ec {Nl+N2)'^ + ^; X)/ -^/(-^z — 1). (for 
simplicity, we assumed an orbital independent Hubbard-like repulsion U for double occupancy 
of a MO). Aea; is a Hund's rule triplet-singlet splitting for the two electronic levels. In 
comparison to Eq. (^ we have neglected single electron hopping mediated by two-electron 
screening effects (which are small for most tt systems in the chosen basis). 

The bias is dropped symmetrically over the molecule, so no capacitive shifts of energies 
appear (this can be easily included). A likely energetic 'term scheme' for the single- and 
two-particle states described by model Hamiltonian Eq. (||) is depicted in the left panel of 
Fig. |l|. The electrodes are considered as non-interacting electron reservoirs. The reservoirs 
are assumed to be occupied according to an equilibrium Fermi distribution function /^(a;) = 
/(w — /Za), where /Zq denotes the electrochemical potential of electrode a = L,R. The MOs 
couple to the electrodes via tunnelling contacts with (possibly) very different coupling strength 

h ■ 

Hmol-leads = (^)'^' {tT4a('^<^o. + ^-C-) , (3) 
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Fig. 1 - Equilibrium molecule energies for model Hamiltonian and resulting current-voltage char- 
acteristics. We choose energies such that the triplet states (T) lies below the singlets (512 Si S2) 
and between the two possible doublets, ei — —12, 62 — —11, (Ae = 1), U = 2, Ae^x ~ 0.5, Ec ~ 4, 
r = 0.004, T = 0.025. 



where pe is the density of states (assumed constant) of the non-interacting electrons in the 
leads, described by operators ako-a. T denotes the scale of the broadening of the MOs due to 
the coupling to the leads. The dimensionless tunneUing parameters t" can depend on a, I. 

In order to investigate whether the predicted effects are robust, we include a couphng of 
the molecule to a (broad band) boson field which simulates the relaxation of excited states in a 
real molecule by coupling of the electrons to an electromagnetic field (photons) and vibrations 
(phonons) pO| . 

Theoretical Approach. We use a Master equation approach for the occupation proba- 
bihties Pg of the molecular many-body states. The transition rate E^s' from state s' to s 
is computed up to linear order in F using golden rule (second order perturbation theory) 
in both the electrode-molecule and the bosonic coupling. For the transition rates we have 
Ess/ — iJ2a p=± ^sf') + ^L' where E"j^ is the tunnelling rate to/from electrode a for creation 
(p=+) or destruction (p=-) of an electron on the molecule. We have 

i:':^ ^TUEg - Eg,)j2\j2tr < s\ci\s' > p, (4) 

and a corresponding equation for S",^ by replacing fa ^ ^ — fa- The boson- mediated rates 
E^j,, describe absorption and emission of bosons. For photons we have 

4g2 

SL' = 9ph—^iE.s - E,,fm{E, - E,,) I < s\d\s' > |2 , (5) 

where d is the dipole operator and Nh{E) denotes the equilibrium Bose function, gph is a 
parameter that allows us to modify the strength of the coupling to simulate increased dipole 
moment or other sources of relaxation, gph = 1, unless noted. This value corresponds to a 
dipole of charge e and length lA. 

We determine the Ps by solution of the stationarity condition: 



Ps — — ^ ^ (Ess'Ps' ~ "^s' sPs) 
s' 



(6) 
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Fig. 2 - Left Panel: I-V-characteristics for various coupling t^. A pronounced NDC effect is observed 
for reduced t2- Right Panel: Occupation probability Ps of the relevant molecule states for t2 = 0.03. 
The fat solid line indicating reaches nearly unity in the blocking regime at bias l^^jg^g > 6. We 
multiply the probabilities with the corresponding degeneracy, so the sum of the Ps adds up to unity. 



The current in the left and right electrode can then be calculated via 

/.=e^(E^+P,,-S^-p,). (7) 

ss' 

The bosonic transition rates do not contribute directly to the current, since they do not change 
the particle number on the molecule. 

Results. The effective Hamiltonian Eq. (|^) affords several generic scenarios for NDC. The 
NDC is generic in the sense that NDC will occur at some bias for an initially charged molecule 
(case (1)) as well as an initially empty molecule (case (2)). 

Case (1): The right panel of Fig. ^ shows the I-V-characteristics for equal tunnelling 
couplings tf^l with a symmetric bias, fiL = ^fJ-R- There are four characteristic steps which 
are related to the onset of the triplet and the three singlet states. From the plateau widths 
all characteristic energy scales can be deduced. 

Strong NDC behaviour is observed if one MO couples much more weakly to the right side 
than the other MO, e.g. = 0.03; = = t2 = 1, see Fig. ||. For a certain bias region the 
current is suppressed by a factor (i|^)^. The reason for this current decrease is the occupation 
of a molecule state (^2 in this case) from which the molecule can not escape anymore due to 
a combination of blocking Fermi sea, Coulomb blockade and the small coupling of an MO to 
the electrode. 

Initially, the molecule is singly occupied in state Di . The current starts at a bias when the 
first two-electron state (triplet) becomes occupied (the "empty" state has higher energy for 
the given parameters). The current can flow via sequential hops through MOi. The electron 
on MO2 is essentially stuck since its tunnelling time to the right reservoir is suppressed by a 
factor (^2^)^, and tunnelling to the left is suppressed because of the blocking Fermi sea. 

But at larger bias the electrons tunnelling onto the molecule from the left can also form 
the state ^2, with both electrons in MO2 as depicted in Fig. ^ No other electron can enter 
the molecule at this bias because of the charging energy. Since the relaxation due to the boson 
coupling is very slow, the only relevant decay of this state is via the small coupling to the 
right electrode. Consequently, the molecule is stuck for a long time in state 5*2. As the right 
panel of Fig. ^ shows the average probability is nearly unity. A relative suppression of 
t2 by 0.3 is sufficient to achieve a pronounced NDC effect. Increasing the temperature will 
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Fig. 3 - Left Panel: I-V-characteristics for various coupling . Again, NDC is observed for reduced 
tf, this time involving as the blocking state. Right Panel: Occupation probability Ps of the relevant 
molecule states for — 0.03. The dashed line indicating Ps^ is nearly unity at bias V|-,jg^g > 4. 



broaden the plateau steps and shift the current maximum slightly to larger bias (not shown). 
At much larger bias (not shown), states with an additional electron become occupied, and the 
current rises again. 

The left panel of Fig. |^ shows that for the same set of energy parameters NDC is observed 
also if if is suppressed instead of if"- In this case the blocking state is Si as indicated by the 
occupation probabilities in the right panel of Fig. ^. 

Case (2): NDC is also observed if we start from an initially uncharged molecule, see the 
left panel of Fig. ^. Note that NDC is observed for — 0.03 for positive bias only, whereas for 
negative bias there is a simple step (dashed curve). This is because if the bias is negative MO2 
and consequently D2 will not become occupied at all, so the current through Di continues to 
flow. We also want to point out that if both the left and right coupling of an MO is reduced, 
NDC is observed for both signs of bias, however, it is limited to a Peak to Valley Ratio (PVR) 
of 2 1 21 1, whereas the PVR for the one sided suppression is limited by 1/(^2')^ ^ 1- 

In the right panel of Fig. ^ we show the influence of internal molecule relaxation by 
increasing the boson coupling gph for case (2) (case (1) shows similar behaviour). An increase 
in the bosonic relaxation rate by six orders of magnitude over the one obtained in dipole 
approximation is necessary to completely eliminate NDC behaviour. It is debatable whether 
coupling to vibrations of the molecule can provide such a rate. However, even in a situation 
where the coupling gph is nominally large there can be selection rules that prevent decay of 
certain states. An example would be the inhibition of (direct) transitions between states of 
different total spin, i.e singlet-triplet transitions. Then, NDC will take place if the triplet 
states (r) are lower in energy than the singlets {Si, S12, S2) and Si is the singlet of lowest 
energy, different from the situation shown in Fig. |l[ The blocking state is Si, the singlet of 
lowest energy. Because of the energy balance "decay" (emission of a boson) of this state could 
only involve the triplets. But this is "forbidden" by the different total spin of the triplets. 
Therefore, as long as there is no absorption of bosons to speak of, this scenario of NDC is 
stable even in the presence of a boson coupling. The rate of absorption processes will strongly 
depend on the number of bosons present. The number of bosons increases with temperature, 
therefore the NDC effect will in general decrease as the temperature is increased due to fact 
that the molecule can escape the blocking state by absorption of bosons. 

The proposed scenarios for NDC are possible explanations for a recent experiment 
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Fig. 4 - I-V-characteristics of the initially neutral molecule, ei — —0.5, £2 ~ 0.5, U = 1.5, Aex = 0.5 
and Ec = 1.5. Left Panel: NDC is observed for t2 = 0.03 involving D2 as the blocking state. If 
both t2 and ^2 ^-r^ suppressed the size of the NDC is limited (see text). Right Panel: Relaxation 
by photon emission gph destroys NDC, but only if the coupling Qph is increased by several orders of 
magnitude over the dipole approximation. 



showing a PVR of more than 1000 at T ~ lOOA'. But in that experiment one is probing 
transport through a whole array of molecules and the nature of one contact side is unknown. 
The temperature dependence is strong (PVR ^ 1.5 at room temperature) and shows a peak 
shift to smaller bias with increasing temperature. Whereas the decrease of the PVR could be 
understood by the increased absorption of bosons (see discussion above) the shift is without 
the reach of the considered model. 

Physical Realization. The energetic arrangement of two closely spaced molecular orbitals 
that couple very differently to the left and right electrodes can be realized in principle with 
very simple aromatic molecules that are suitably substituted to break their full symmetry. As 
an example. Fig. ^ shows the LUMO and LUMO+1 for a 1,3-dimethyl benzene (meta- xylene) 
that are closely spaced energetically and that couple very differently on various possible un- 
substituted contact sites. There will be MOs which are antisymmetric with respect to the 
2,5-axis mirror symmetry with vanishing wave function amplitude at the 2 and the 5 position. 
In contrast, the symmetric MOs will in general have non- vanishing wave function at these po- 




Fig. 5 - The LUMO and LUMO+1 for a double methyl substituted benzene. By (anti)-symmetry, 
one of the MOs will have no coupling at the 2 position. 
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sitions. If one couples the molecule at the 2 and 6 positions to electrodes, the LUMO couples 
to both electrodes whereas the LUMO+1 would have no coupling to the electrode 'connected' 
to the 2 position |^^. Thus, the situation of a strongly MO and electrode dependent coupling 
seems to be generic for small aromatic molecules with ligand groups. 

Conclusions. We have developed a model of non-linear charge transport through a metal- 
molecule-metal nanojunction. We have shown that the interplay of charging effects and 
the spatially non-trivial electronic structure of the molecule can lead to current peaks and 
strong negative differential conductance. NDC will be observed no matter whether the tunnel 
coupling of the lower or the higher MO is suppressed. For a coupling to photons in dipole 
approximation the relaxation rate induced by the photons is several orders of magnitude too 
small in comparison to typical tunnelling rates to have an effect. We believe that the model is 
sufficiently generic to be realized in certain classes of aromatic molecules with tunnel contacts 
to electrodes. 
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